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ABSTRACT 

The  formulation  and  solution  of  inverse  problems  for  the  estimation  of 
parameters  which  describe  damping  and  other  dynamic  properties  in  distributed 
models  for  the  vibration  of  flexible  structures  is  considered.  Motivated  by  a 
slewing  beam  experiment,  the  identification  of  a  nonlinear  velocity  dependent 
term  which  models  air  drag  damping  in  the  Euler-Bernoulli  equation  is  investi¬ 
gated.  Galerkin  techniques  are  used  to  generate  finite  dimensional  approxima¬ 
tions.  Convergence  estimates  and  numerical  results  are  given.  The  modeling 
of,  and  related  inverse  problems  for  the  dynamics  of  a  high  pressure  hose  line 
feeding  a  gas  thruster  actuator  at  the  tip  of  a  cantilevered  beam  are  then 
considered.  Approximation  and  convergence  are  discussed  and  numerical  results 
involving  experimental  data  are  presented. 
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L.  Introduction 


The  purpose  of  this  note  is  to  illustrate  and  explain  some  of  the  ideas  underlying  the  use  of 
parameter  estimation  techinques  in  investigating  damping  and  other  dynamic  phenomena  in  several 
classes  of  distributed  models  for  flexible  structures.  We  do  this  in  the  context  of  two  specific 
examples  drawn  from  experimental  structures.  In  the  first  example  we  present  techniques  and  results 
that  can  be  used  to  study  nonlinear  aspects  of  viscous  damping  (nonlinear  air  drag)  in  slewing 
maneuvers  with  flexible  beam  like  structures.  We  shall  describe  some  fundamental  questions, 
present  a  model  for  which  an  estimation  problem  is  of  importance,  and  then  show  how  this  inverse 
problem  can  be  approximated  for  computational  purposes.  We  do  this  in  a  weak  or  variational  setting 
and  give  convergence  arguments  to  provide  a  theoretical  foundation  for  the  numerical  schemes  we 
have  used.  These  convergence  results  are  then  followed  by  presentation  of  a  numerical  test  example. 
Our  methods  have  proved  useful  in  current  studies  with  experimental  data,  the  results  for  which  will 
be  presented  in  detail  elsewhere. 

In  a  second  example,  we  outline  techniques  that  have  been  useful  in  developing  accurate  models 
for  the  dynamic  effects  of  a  flexible  gas  hose/tip  mass/thruster  apparatus  when  it  is  attached  to  a 
flexible  beam  to  provide  an  active  control  system  (the  so-called  "RPL  experiment").  Since  detailed 
mathematical  arguments  for  this  project  are  given  elsewhere  [BGRW],  we  shall  only  outline  the 
model,  the  approximation  ideas,  and  present  a  summary  of  numerical  results  obtained  when  using  our 
general  approach  with  experimental  data. 

2L  Nonlinear  Damping  in  Slewing  Maneuvers 

In  a  series  of  papers  [HJ],  [JH],  [JHR],  Juang  and  his  co-workers  describe  experiments  carried 
out  to  demonstrate  the  feasibility  of  actively  controlling  (stabilizing)  vibrations  of  a  beam  during 
slewing  maneuvers.  Experiments  were  conducted  with  a  1  meter  steel  beam  and  with  a  3.9  meter 
aluminum  honeycomb  solar  panel  cantilevered  in  a  vertical  plane  and  rotated  in  the  horizontal  plane. 
Each  was  attached  to  a  torquing  motor  at  the  hub  of  rotation  or  "root"  of  the  beam.  In  typical 
experiments,  the  beams  were  slewed  30°  to  45°  in  1.5  to  4.5  seconds.  Strain  gauges  were  located  at 
the  root  and  at  .22  and  .5  of  the  length  l  of  the  beam.  An  angle  potentiometer  (at  the  root)  measured 
angular  displacements  during  the  slew.  These  measurements  (strain,  yxx,  and  angular  displacement, 
0)  were  used  as  feedback  to  the  motor  which  was  then  used  to  suppress  the  vibrations  via  an  LQR 
theoretical  formulation  for  the  feedback  control  laws. 

A  series  of  slewing  maneuver  experiments  were  carried  out  in  a  laboratory  at  NASA  Langley 
Research  Center  with  possible  effects  due  to  air  damping  present.  These  were  then  followed  by 
repetition  of  the  experiments  in  a  vacuum  chamber.  Experiment  and  theory  were  in  good  agreement 
in  the  vacuum  chamber  experiments,  but  there  were  significant  discrepancies  between  the  theoretical 
model  based  simulations  and  the  experimental  data  obtained  in  the  laboratory  setting.  It  is  important 
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to  understand  the  variation  in  responses  of  controlled  flexible  structures  in  such  differing 
environments  since  many  future  large  spacecraft  studies  will,  by  necessity,  involve  model 
extrapolation  and  design  based  on  laboratory  performance  only  without  the  benefit  of  vacuum 
chamber  comparisons. 

It  has  been  suggested  that  the  inaccuracies  in  the  Juang,  et.  al.  simulations  most  likely  resulted 
from  two  types  of  model  error:  (1)  the  absense  of  a  viscous  damping  term  to  represent  the  nonlinear 
air  damping  and  (2)  dissipation  at  the  fixed  end  of  the  cantilevered  beams  was  not  included  in  the 
model  even  though  there  was  some  loose  "play"  in  the  clamped  or  "built-in"  end  of  the  beams.  This 
latter  mechanism  for  energy  absorption  presumably  should  be  modeled  by  some  type  of  nonlinear 
boundary  condition  in  place  of  the  usual  no  displacement,  no  slope  boundary  conditions:  y  =  yx  =  0. 


Figure  2.1 

Here  we  focus  on  a  nonlinear  air  damping  component  of  the  model  as  suggested  by  Juang,  et. 
al.,  and  show  that  one  might  effectively  compute  this  using  the  experimental  data  in  a  least  squares 
setting.  We  use  a  formulation  proposed  by  Juang  and  his  co-workers;  this  model  can  be  derived  by 
combining  first  principle  energy  considerations  with  laboratory  findings  on  nonlinear  drag  forces.  If 
0(t)  and  y(t,x)  denote  respectively  angular  displacement  (from  some  reference  angle)  at  time  t  and 
bending  deflection  along  the  beam  at  time  t  and  position  x  as  depicted  in  Figure  2.1  above,  the  model 
including  actuator  dynamics  has  the  form  (see  [JHR]  for  a  more  detailed  discussion) 


•V. 

tv\ 

•V', 


d2y  d2e  de  dy  ,  dd  dy ,  d\ 

at)  p-r  +  px—  *  C3f(x_  +  X)  4.  C„(X_  4  x)  +  El  _  =  0 


_  ,  d  0  r*  dy  ,  ^  d0  dy.  ,  r<  d9  9y.  , 

(2.2)  In -  +  f  px — -  dx  +  I  c,f(x—  +  -~)  xdx  +  c.{x—  +  xdx  =  T(t) 

B  dt2  J°  ^2  Jo  3  dt  9t  ■’o  4  dt  9t 
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Here  et  represents  the  applied  armature  voltage  (this  term  involves  strain  feedback  in  the  control 
problem),  the  terms  involving  c3  and  c4  represent  viscous  damping  and  the  coefficients  p  and  El  are 
the  usual  beam  material  parameters,  linear  mass  density  and  flexural  stiffness,  and  IB  is  the  moment 
px2dx  about  the  axis  of  rotation.  The  damping  function  f  is  assumed  to  have  the  form  f(v)  =  vlvl 
for  v  in  some  bounded  region.  Of  course,  one  can  assume  without  loss  of  generality  based  upon 
physical  considerations  that  f  becomes  bounded  as  v  becomes  large. 

A  number  of  state  space  formulations  of  this  model  can  be  investigated  in  the  context  of 
identification  and  control  problems  including  the  following: 

•  d0 

(i)  Both  0  and  0  =  —  are  treated  as  states, 

dt 

(ii)  The  time  history  for  0  is  assumed  known,  but  0  is  treated  as  a  state  (if  0(t)  is  known 
from  experimental  data,  it  does  noi  follow  that  0(t)  can  be  effectively  obtained  by 
differentiation  of  the  data), 

(iii)  Both  0  and  0  are  measured  reliably  so  that  each  can  be  treated  as  a  known  quantity. 

For  our  discussions  here  of  estimation  of  the  damping  coefficients,  we  assume  that  (iii)  holds  so 
that  one  can  combine  (2.1)  and  (2.2)  to  also  eliminate  0  as  an  unknown  in  the  model.  We  do  this 
before  formulating  a  least  squares  problem  for  estimation  of  the  damping.  We  also  include  a 
Kelvin- Voigt  material  damping  term  in  the  beam  dynamics  equation. 

We  consider  then  the  coupled  system  (actuator  dynamics  plus  transverse  vibrations  of  the  beam) 
for  0<  x  <  l ,  t  >  0 

(2.5)  pytt  +  px0  +  c3f(x0  +  yt)  +  c4{x0  +  yt)  +  El  D4y  +  cDID4yt  =  0, 

(2.6)  IB  0  +  pxyn  dx  +  c3  f(x0  +  yt)  xdx  +  c4  J  {x0  +  yt}xdx  =  k0  eg  -  kj0  -  k  20, 
with  appropriate  boundary  conditions  for  1 1 0 

(2.7)  y(t,0)  =  Dy(t,0)  =  0, 

(2.8)  EIDVt,*)  +  cDID2y,(t,i)  =  0,  El  D3y(U)  +  cDI  D^t,*)  =  0, 

and  appropriate  initial  conditions.  Here  and  below  we  shall  use  the  notation  D  =  9/9x.  Assuming 
that  0(t)  and  0(t)  are  known  for  0  <  t  <;  T,  where  T  is  the  duration  of  the  experiment,  we  may  solve 
for  0  in  equation  (2.6),  subsitute  this  into  equation  (2.5),  and  obtain  a  single  nonlinear  partial 
differential  equation  for  the  bending  deflection  y.  Upon  doing  this  we  obtain 


(2.9) 


where 


pytt  +  U.ytt)  +  c3  7X  (t,x,yt)  +  c4  72  (t,x,yt)  +  g(t,x) 

+  c3f(x0  +  yt)  +c4{x0  +  yt}  +  El  D4y  +  cDI  D4yt  =  0 


•VfiY ivtY*'  *;  vv  v  V. 


I 

'  f 


V 

tv4 

,vJ 

'It 


(2.10) 

& 

III 

1 

(2.11) 

5j(t,x,yt) 

=-  iBp+k2i„f<se+>,.>sds' 

(2.12) 

J2(t,x,yt) 

= •  iBp+k2i„|s9+y.|sds- 

and 

(2.13) 

g(t,x)=  - 

1-^.(k0e,<o-k1e(t)). 

The  boundary  conditions  (2.7),  (2.8)  involving  no  displacement,  no  slope  at  x  =  0  and  no  moment, 
no  shear  at  x  =  i ,  can  be  equivalently  written  as 

(2.14)  y(t,0)  =  0,  Dy(t,0)  =  0,  DtyU)  =  0,  D3y(tJ)  =  0. 

Using  this  system,  we  may  formulate  an  estimation  problem  for  the  parameter  q=  (qpq2.Q3.q4) 
=  (El,  cDI,  c3,c4),  given  observations  e;,  i  =  1,  ...,  M,  for  the  root  strain  D2y(ti,0).  The  problem 
becomes  one  of  choosing  from  some  admissible  set  of  parameters  QcR4  (we  treat  here  only  the 
constant  parameter  case  -  variable  coefficients  can  be  readily  treated  with  appropriate  modifications  of 
the  arguments  -  see  [BC1],  [BC2],  [BCR],  [BR1]  for  details)  a  parameter  q*  which  minimizes  over 
Q  the  least  squares  criterion 

M 

J(q)=  X 1  Dy^.o-.q)  -  1 

i=l 

where  y(-,  • ;  q)  is  the  solution  of  (2.9),  (2.14)  for  a  given  set  of  initial  conditions  y(0,x)  =  d>(x), 
yt(0,x)  =  'F(x).  Least  squares  problems  of  this  type  possess  a  number  of  interesting  aspects  (infinite 
dimensional  states,  theoretical  and  computational  ill-posedness,  etc.)  which  have  been  discussed  by 
us  and  others  in  a  number  of  previous  publications.  Here  we  focus  on  one  particular  question:  state 
approximation  techniques  and  the  convergence  arguments  for  approximating  parameter  estimates 
(these  convergence  arguments  are  also  an  important  part  of  theoretical  results  which  guarantee  a  type 
of  inverse  method  stability  -  i.e.  continuous  dependence  of  estimates  on  the  observations  -  see  [B]). 
Before  turning  to  these  convergence  arguments,  we  rewrite  the  system  in  a  variational  form 
(conservative  form)  with  states  v  =  yt,  w  =  D2y. 


We  first  rewrite  (2.9)  as 


(2.15) 

(2.16) 
(2.17) 


pv,  +  Z(vt)  +  q3  7,0, x,v)  +  ^(t.x.v)  +  g 

+  q3f(x8  +  v)  +  q4[x0  +  v)  +  q,  D2w  +  q2  D4v 


=  0 


w{  =  D\ 
v(0)  =  4/, 


w(0)  =  D2^, 
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where  we  seek  solutions  v  e  H4  n  ,w  €  H|  with 

H2  =  {<pG  H2(0,£)  I  (p(0)  =  D9(0)  =  0} 

Hr  =  {9  e  H2(0,-6)  I  9(^)  =  D9(4)  =  0). 

Then  an  equivalent  weak  or  variational  form  is:  Find  (v,w)  e  H[xH°  satisfying  the  initial  condtions 
(2.17)  and 

<2-18)  <  pvt  +  L(vt)  +  q3J1(t,-,v)  +  q4J2(t,-,v),  9  > 

+  <  g(t.O  +  q3f(x0  +  v)  +  q^x©  +  v),  9  >  +  <  q}w  +  q2D2v,  D29  >  =  0  for  9  e 

(2.19)  <  wt- D2v,\)/>  =  0  for  ye  H°. 

We  next  make  the  following  observations:  for  each  of  the  terms  involving  q*  in  equation  (2.18), 
there  is  an  analogue  involving  q3,  the  difference  being  that  the  q3  terms  are  nonlinear,  the  q4  terms  are 
linear.  Since  our  primary  emphasis  here  is  in  presenting  convergence  arguments  for  estimation  of  the 
nonlinear  damping  coefficients  and  since  these  arguments  readily  extend  to  the  case  involving  the 
linear  q^  terms,  we,  for  the  sake  of  exposition,  shall  assume  throughout  the  remainder  of  our 
discussions  that  q4  =  0.  Thus  the  system  of  interest  to  us  is 

(2.20)  <  pvt  +  £(vt)  +  q3^1(t,,,v)  +  q3f(x0(t)  +  v)  +  g(t,  ),  9  > 

+  <  qjw  +  qjD^,  D29  >  =  0  for  9  e  , 

(2.21)  <wt-D2v,\jr>  =  0  foryeH0, 

(2.22)  v(0)  =¥,  w(0)  =  D20>. 

We  approximate  the  system  via  the  usual  Galerkin  procedures  by  choosing  finite  dimensional 
subspaces  c  H£  and  HN  c  H°  satisfying  certain  approximation  properties  (as  N  -»  «>)  to  be 
specified  below.  The  approximate  system  for  vN  e  H^,  wN  e  HN  is  given  by 

(2.23)  <  pv“+  )  +  q3J1(t,-,vN)  +q3f(x0  +  vN)  +  g(t,),  9N  > 

+  <  qjWN  D29n)  =0  for  9N  e  HJJ , 

(2.24)  <  w^-  D^14,  \j/N  >  =  0  for  yN  e  HN, 

(2.25)  vN(0)  =  P£v(0),  wN(0)  =  PNw(0) 

where  P£,  PN  are  the  orthogonal  projections  (in  the  H°  norm)  of  H°  onto  H£,  HN,  respectively.  We 
then  may  define  a  sequence  of  approximating  parameter  estimation  problems  consisting  of  minimizing 
over  Q  the  criterion 

(2.26)  JN(q)  =  X 1  wN(t,0;q)  -  e.l2 

i=l 

where  wN  is  the  solution  of  (2.23)  -  (2.25)  corresponding  to  q. 


The  procedure  for  arguing  convergence  of  the  approximate  parameters  qN  (a  minimizer  of  JN  in 

(2.26) )  is  now  well-documented  in  a  number  of  our  previous  papers  (see  [B]  for  a  summary).  A 
crucial  step  entails  the  following:  If  {qN}  is  an  arbitrary  sequence  in  Q  with  qN  — >  q  in  Q,  one  must 
establish  the  convergence  zN  -» 0,  uN  — >  0  where  zN(t)  =  vN  -  P£  v,  uN(t)  =  wN  -  Pw  with  vN,  wN  the 
solutions  of  (2.23)-(2.25)  corresponding  to  qN  and  v,w  the  solutions  of  (2.20)  -  (2.22) 
corresponding  to  q.  We  outline  the  arguments  to  establish  this  convergence. 

Using  (2.20)  -  (2.22)  and  (2.23)  -  (2.25)  with  q  =  qN,  we  find  for  all  9  e  H£ 

<pz?,9>  =  <vN-vt  +  (I-P^)vl)(p> 

=  <  (I  -  P£ )  vt  +  £(vt)  -  £{vN),  9  > 

+  <  q3^(t,-,  v)  -  q^O,  •,  vN)  +  q3f(x0  +  v)  -  q^f(x©  +  vN),  (p  > 

+  <  qjW  -  q^wN  +  q2D2v  -  q^D^14,  D2  (p  > 

=  <  (I  -  P£)vt+  K-z*4)  +  £((I  -  P{J  )vt),  9  >  +  (q3  -  q3  )  <  71(t,-,v),  cp  > 

+  q3  <  ^(t,-,v)  -  J^t,^),  cp  >  +  (q3  -  q^)  <  f(x0  +  v),  (p  > 

+  q3<  f(x9  +  v)  -  f(x0  +  vN),  <p  >  -  <  q^(wN  -  Pw )  +  q£D2(vN  -  P£  v),  D2  cp  > 
+  <  qtw  -  q^PNw  +  -  q2  D2  P£  v,  D2  <p  >  . 

Thus  we  find  (suppressing  some  of  the  obvious  function  arguments) 

<pzN+£(zN))(p>  =  <(I-P£)vt  +  U  (I  -  P£  )v,),  cp  > 

+  (q3  -  q")  <  TX(V),  cp  >  +  <  fx(v)  -  7j(  P£  v),  <p  > 

(2.27)  +  q3  <  P£  v)  -  ?\(vN),  cp  >  +  (q3  -  q^1 )  f(x0  +  v),  cp  > 

+  q^  <  f(x0  +  v)  -  f(x9  +  P£  v),  cp  >  +  q^  <  f(x0  +  P£  v)  -  f(x0  +  vN),  cp  > 

-  <  q^uN  +  q^D2zN,  D2  9  >  +  <  qjW  -  q^  PNw  +  qjD^  -  q2  D2  P£  v,  D29  > 

for  all  9  e  H£.  Also,  for  y  e  HN  we  find 

<  u*4,  Y  >  =  <  wtN  -  PNwt,  y  >  =  <  w*4  -  wt  +  (I  -  PN)  wt,  y  > 

(2.28)  =  <  -  D^,  y  >  +  <(I  -  PN)wt,  Y  > 

=  <  D2zn,  y  >  +  <  D2(  Pl  -  I)v,  Y  >  +  <  (I  -  PN)w,,  y  >• 

To  obtain  convergence  estimates  we  make  particular  choices  for  9  and  y  <n  (2.27)  and  (2.28).  Let 
9  =  zN  e  H£  in  (2.27)  and  y  =  uN  e  HN  in  (2.28)  and  define  0*  h  D2(  P£  -  I)v,  0£  s  (I  -  PN)wt. 
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Then  from  (2.28)  we  obtain,  using  the  inequality  ab  <  —  a  +  eb  where  e  >  0  can  be  arbitrarily 
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chosen. 


(2.29) 
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2  dt  4e  1  4e  2  4e 
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=  e  l  D  z  l  + —  I  u  I  +  e  10.  I  +  e  10,  I  . 

4e 

Before  deriving  an  analogous  estimate  from  (2.27),  we  require  some  hypotheses  on  the  nonlinear 
damping  function  f  which  will  in  turn  yield  estimates  on  L  and  7V  Recall  that  we  can  expect  f(v)  to 
behave  like  vlvl  =  (sgn  v)  v2  for  /  in  some  neighborhood  (not  necessarily  small)  of  the  origin,  while  it 
becomes  bounded  for  v  large.  Thus,  for  our  theoretical  considerations  here  we  make  the  reasonable 
assumptions  that  f  satisfies  a  Lipschitz  condition 


(2.30) 


lf(^)  -  f($2)  I  *  KIV^I 


for  real 


as  well  as  satisfying  a  boundedness  condition 


(2.31) 


I  f(£)  I  <  M. 


From  (2.30)  we  readily  obtain  the  estimates  (we  shall  need  these  to  use  in  estimates  from  (2.27) 
with  (p  =  zN) 


(2.32) 


<  f(x0  +  v)  *  f(x0  +  P^  v),  zN>  ^  JQ  K  1  (I  -  P^)  v(t,x)  1 1  zN(t,x)  I  dx 
<  K  I  (I  -  P^  )v(t)  1 1  zN(t)  I  <  |  I  (I  -  P*v(t)  I2  + 1  I  zN(t)  I2, 


where  we  have  used  I  I  to  denote  either  the  absolute  value  or  the  H°(0,£)  norm  (the  interpretation 
being  clear  from  the  usage).  Similarly,  we  find 


(2.33)  <  f(x0  +  P£v)  -  f(x0  +  vN),  zN  >  ^  K  I  vN  -  P[Jv  1 1  zN  I  =  K  I  zN(t)  I2. 


Furthemore,  using  the  definition  (2.1 1)  of  7V  we  easily  find 


(2.34) 


J,(t,-,v)  -  J  (t.-.P^v),  zN  >  <  J  ~X  f  K  I  (I  -  P^)v(t,s)  I  s  d  s  I  zN(t,x) 
1  i  l  -»o  IB  +  kj  Jo  L 

<  p[l(I-P|Vl2  +  lzN(t)l2], 


where  the  constant  p  depends  in  an  obvious  way  only  on  l  ,K,p,  IB,  k2.  With  similar  calculations 
we  find 


(2.35)  <  7,(t.-,  P£  v  )  -  7,(t,  •  vN),  zN  >  <  2p  I  zN(t)  I2. 


Using  the  definition  (2. 10)  of  L  we  next  find 
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N.  N  f* 
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(2.36) 


^/‘ps.MdsjzV 

_£L  I  iff‘xZN(.,x)dxl 

iB  +  kj  2  d.  [Jo  J 


-p  Id  N.  .  2 

T  —  <  X,  Z  (t)  >  . 


Ig  +  l^  2  dt 


Also 


(2.37) 


<  a  (I  -  p“)  vt),  ZN  >  <  jM—  I  (I  -  p")  v(  1 1  z"  I 

B  z 

<  U  [  I  (I  -  Pl  )  vt  I2  +  I  zN(t)  I2] 
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where  the  constant  \x  depends  only  on  p,  IB,  k^,  l . 

Finally,  we  can  now  obtain  a  desired  estimate  from  (2.27)  with  cp  =  zN.  Using  (2.32)  -  (2.37) 
we  find 


1  £ 

2  dt 


.  N  .2  p  N  2 

p  I  Z  I  -  - - f-f—  <  X,  Z  > 


*B  +  ^2 

+  ii[  KI-PlK  l2  +  lzN|2]  +  |q3.qJ|  M  I  zN  I 
(2.38)  +  I  q"  I  p[  I  (I  -  )v  I2  +  I  zN  l2J  +  |q*  |  2p  I  zN  I2  +  lq3  -  q*  I M JT  I  zN 

+  lq3  I  j  [  I  (I '  )v  I2  +  I  zN l2j  +  |<£|  K  tzN|2 

+  1  q*  1 1  uN  1  !D2zN  l  -  <  q"  D2zN,  D2ZN  >  +  10^1 1D2ZNI 

where  8^  =  (qtw  -  q^PNw  +  q2D2v  -  q^D2P£  v). 

N  /  3  m 

Noting  that  <  x,  z  >  <  yj  t  / 3  I  z  I,  we  find  that 


(2.39) 


2  2  3 

-P  N  2  -p  i  ,  N.2 

— —  <x,z  >  >  — - —  — 


TB  +  k2 


IB+k2  3 


i  z  r  = 


-p 


N  ,2 


1  + Ikjpl 

ft  2  3 

since  I  =  px  dx  =  p  £  n  We  assume  that  Q  is  a  compact  set  and  that  the  admissible  parameter 
set  Q  entails  the  constraints  q2  >  v  for  some  positive  constant  v.  We  then  have 
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v  I  D2zn  I2  <  <  q^  D2zn,  D2zn  >. 


Using  these  estimates  in  (2.38)  we  obtain 


Id  -  N  .2  n2  N  2  1  M  - 
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or,  assuming  that  I  q3  -  q^l  — »  0, 

(2.40)  T  ^  P  1  zN  |2  -  <  x>  zN  >2  <  3^1  ZN  |2  +  ^ (t)  +  (2e  -  v)  I  D2zN  I2 

.2 

,  k  N  .2 
+  —  I  u  I 

4e 

wnere  /?*  (t)  is  bounded  and  -»  0  as  N  — »  °o  under  the  usual  assumptions  (see  [BJ,  [BCR],[BR1])  on 
the  approximation  properties  of  and  HN  (i.e.  that  P£  -»  I,  PN  ->  I  in  the  desired  topologies).  Foi 
example,  approximations  based  upon  cubic  splines  in  H£  and  linear  splines  in  H°  would  suffice. 
Finally  combining  (2.29)  and  (2.40)  we  obtain 


Id  -  N  ,2 

(2.41) 


h  +  *2 


N  2  .  N.2 

<  X,  z  >  +  I  U  I 


£  (3e  -  v)  I D  z  I 


+  Yj  I  zN  I2  +  Y2luNl2  +  GN(t) 

where  GN  is  bounded  with  GN(t)  — »  0.  Integrating  this  inequality  and  using  the  fact  that  zN(0)  =  0, 
uN(0)  =  0,  we  find 

I  p  I  zN(t)  I2  -  j-fi-  —  <  x,  zN(t)  >  2  +  I  uN(t)  I2  +  (V .  3e)  |  D2zN(s)  I2  ds 

(2.42) 

<  JQ  [yi  I  zN(s)  I2  +  Y2  I  uN(s)  l2J  ds  +  J^[lGN(t)  I  +  I  At)  l]  dt . 


Finally,  using  (2.39)  and  choosing  e  such  that  v  -  3e  =  8  >  0  in  this  last  inequality,  we  obtain 
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1  +  3k2/pe' 

+  8  I  D2zN(s)  I2  ds  <  A(N)  +  J*  [yi  1  zN(s)  I2  +  Y2  •  As)  l2j  ds 


where  A(N)  — »  0.  Thus,  by  the  usual  Gronwall  arguments  we  have  I  zN(t)  I2  — >  0,  I  uN(t)  I2  — »  0, 
and  ft  I  D2zN(s)  I2  ds  — »  0,  where  zN(t)  =  vN(t)  -  v(t)  and  uN(t)  =  wN(t)  -  PNw(t).  Under 
appropriate  convergence  properties  for  P£  and  P1^,  we  may  then  use  the  triangle  inequality  to  obtain 


I  vN(t)  -  v(t)  I  — >  0, 1  wN(t)  -  w(t)  I  ->  0,  f  I  D2vN(s)  -  D2v(s)  I2  ds  ->  0. 

Jo 


We  note  that  the  above  arguments  can  be  used  to  establish  convergence  at  each  t  of  the  strain  in 
the  H°(0,-fc)  norm.  If  we  use  the  root  strain  D2y(t,0)  in  the  least  squares  criterion,  a  stronger  strain 
convergence  (pointwise  in  the  spatial  coordinate)  is  needed  to  complete  the  theoretical  development. 
Arguments  in  the  spirit  of  those  given  above  can  be  made  to  give  strain  convergence  in  the  H1  norm 
which,  of  course,  yields  root  strain  convergence.  Arguments  of  this  nature  have  been  given  for 
similar  problems  elsewhere  [BCK],  [BR2];  they  involve  some  technical  detail  and  we  shall  not 
pursue  the  development  here.  Instead  we  turn  to  a  brief  discussion  of  some  computational  aspects  of 
these  schemes. 

The  approximation  and  estimation  schemes  discussed  above  can  be  readily  used  to  develop 
computational  algorithms  for  estimation  of  damping  coefficients  (including  those  for  the  nonlinear 
viscous  damping  terms).  We  have  developed  and  tested  numerically  some  software  packages  based 
on  these  ideas.  We  are  currently  using  the  packages  with  experimental  data  provided  by  J.  Juang; 
these  results  will  be  reported  elsewhere.  We  close  this  section  with  a  brief  summary  of  findings  for 
one  of  the  numerical  test  examples  we  have  investigated. 

A  test  example  was  considered  for  the  system  (2.9)  with  l  =  p  =  El  =  1.0  assumed  known  and 
c4  =  cD  =  0.  We  sought  to  estimate  the  damping  coefficient  c3  from  simulated  data  for  the  strain. 
That  is,  we  chose  a  particular  function  y(t,x)  =  .5t2(x4  -  4x3  +  6x2)  as  the  true  solution  of  the  system 
equation  (2.9)  with  an  appropriate  forcing  function  added  to  the  equation.  The  true  parameter  value 
c3  =  2.0  was  used  in  (2.9)  along  with  choices  of  0(t)  and  ea(t)  which  were  qualitatively  similar  to  the 
corresponding  experimental  time  histories  available  to  us.  For  example,  we  used 


9(t)  =  41 


1  -cpy 

1 - e  sin(wn  J3t  +  <j>) 


with  P 


4 


1-C2, 


4>  =  Arctan  P/£ 


and  £,  con  chosen  so  that  0  has  a  maximum  amplitude  of  4.5  at  t  =  1.75.  The  values  kg  =  52.136,  k, 
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=  1672.96  and  Ig  +  k2  =  32.95  were  used  in  (2.9)  along  with  0(t)  obtained  by  differentiating  0. 

For  the  approximating  system  we  used  modified  cubic  splines  for  both  and  HN;  the  basis  for 
was  the  usual  cubic  B-spline  basis  modified  to  satisfy  the  boundary  conditions  of  while  the 
usual  B-splines  modified  to  satisfy  the  boundary  conditions  of  H|  were  used  to  generate  HN.  In  each 
case  the  subspaces  H£,  HN  had  dimension  N  +  1.  In  the  fit  criterion  we  used  the  root  strain  yxx(tj,0) 
at  29  equally  spaced  observations  in  the  time  interval  0  <  t  <  7.  For  the  initial  guess  c^  =  2.5,  a 
residual  of  1789  is  obtained.  For  N  =  4,  the  scheme  produced  a  converged  estimate  03=1 .996  with 
a  corresponding  residual  of  .0088.  A  number  of  other  test  examples  were  studied  with  equally 
satisfactory  findings. 


The  RPL  structure  is  an  experimental  apparatus  which  was  designed  and  constructed  at  the 
Charles  Stark  Draper  Laboratory  in  Cambridge,  Massachusetts  with  funding  supplied  by  the  United 
States  Air  Force  Rocket  Propulsion  Laboratory  (RPL).  Its  primary  function  is  to  serve  as  a  test  bed 
for  the  purpose  of  investigating  control  algorithms  and  instrumentation  (sensors,  actuators, 
processors,  etc.)  for  the  large  angle  slewing  of  spacecraft  with  flexible  appendages.  It  was  designed 
to  specifically  incorporate  those  features  which  make  control  design  for  large  flexible  spacecraft  an 
especially  difficult  and  challenging  problem.  In  particular,  this  includes  light  damping,  high 
flexibility,  a  large  number  of,  and  closely  spaced  natural  modes  of  vibration,  difficult  to  model  and 
coupled  structural  and  actuator  dynamics  and  dissipation  mechanisms,  etc. 
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The  structure  itself  consists  of  four  aluminum  beams,  each  of  length  l  =  4  feet,  width  b=6 
inches  and  thickness  h  =  .125  inches,  cantilevered  in  a  symmetrical  fashion  to  a  central  hub  which  is 
mounted  on  an  air  bearing  table.  The  air  bearing  table  allows  for  the  near  frictionless  rotation  of  the 
entire  structure  about  the  vertical  axis.  Control  actuation  is  achieved  via  nitrogen  cold  gas  thrusters 
mounted  at  the  tips  of  two  opposing  appendages.  The  other  two  beams  are  passive  with  masses  at 
their  tips  serving  only  to  maintain  the  over  all  symmetry  of  the  structure.  Nitrogen  gas  is  supplied  to 
the  thrusters  from  storage  tanks  mounted  to  the  central  hub  through  stainless  steel,  wire  mesh 
wrapped,  flexible,  high  pressure  hoses.  Electro-mechanical  valves  control  the  expulsion  of  the  gas 
from  the  thruster  nozzles.  Each  appendage  is  instrumented  with  a  linear  accelerometer  at  the  tip.  Data 
from  the  sensors  is  recorded  and  control  input  signals  are  generated  using  a  MINC  11/23 
microcomputer. 

Effective  control  design  depends  heavily  upon  the  availability  of  a  high  fidelity  model  for  the 
plar ..  In  the  case  of  the  RPL  structure,  it  is  immediately  clear  that  a  model  involving  partial 
differential  equations  would  be  of  some  use.  For  the  transverse  vibration  of  the  passive  beams,  a 
distributed  parameter  model  based  upon  the  Euler-Bemoulli  equation  together  with  appropriate 
boundary  conditions  describing  the  coupled  motion  of  the  tip  mass  and  the  rigid  body  rotation  of  the 
central  hub  would  be  adequate.  For  the  active  appendages  (i.e.,  those  with  the  tip  thrusters)  on  the 
other  hand,  a  more  sophisticated  model  which  also  captures  the  coupled  dynamic  effects  (i.e. 
additional  mass,  stiffness  and  dissipation,  torsional  motion,  etc.)  due  to  the  motion  of  the  flexible 
thruster  hoses  is  needed. 

Since  the  transverse  vibration  of  each  of  the  individual  appendages  is  decoupled,  for  our 
investigation  here,  we  consider  the  problem  of  modeling  the  hose  effects  on  the  transverse  vibration 
of  a  single  cantilevered  (i.e.  clamped  -  free)  beam  (see  Figure  3.2). 


Figure  3.2 

We  describe  a  model  which  was  suggested  by  S.  Gates  of  the  Control  and  Right  Dynamics  division 
of  the  Draper  Laboratory  wherein  the  hose  is  treated  as  a  damped  linear  harmonic  oscillator  that  is 
rigidly  attached  to  the  thruster  assembly  at  the  free  end  of  the  beam.  More  precisely,  the  hose  is 
modeled  as  a  proof  mass  which  reacts  against  the  tip  or  thruster  mass  via  an  elastic  spring  and  a 
linear,  viscous  damper  (see  Figure  3.3) 
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Figure  3.3 


Letting  u(t,x)  denote  the  vertical  displacement  of  the  beam  at  time  t  and  position  x,  0  <  x  <  ( , 

and  assuming  only  small  deformations  ( i.e.  Iu(t,x)  I  <  <  (,  I  (t,x)  I  «  1) 

dx 

we  use  the  Euler-Bemoulli  equation  together  with  Kelvin- Voigt  viscoelastic  material  damping  to 
describe  its  transverse  vibration  .  That  is 


-\4 

(3.1)  p  (t,x)  +  cDI  — -  (t,x)  +  El  — ^  (t,x)  =  0  0  <  x  <  (,  t>0 

dt2  dx4  3t  dx4 


where  p  is  the  linear  mass  density,  cD  and  E  are  respectively  the  material  coefficient  of  viscosity  and 
modulus  of  elasticity  and  I  =  bh3/12  is  the  second  moment  or  moment  of  inertia  of  the  uniform, 
rectangular  cross  section  of  the  beam. 

At  the  free  end  of  the  beam,  corresponding  to  x  =  ■(.the  motion  of  the  tip  mass  rrij  and  the 
hose  mass  mH  are  given  by 

d2u  .  d3  du  d3u 

m  (t,()-c  l—  _(t,  ()-EI  — (t,  () 

(32)  dt2  dx3  dt  dx3 

=  CH(^(t)  *  |f(t’  *)}  +  kH(y(t)  '  u(t’*)}  +  f(t)’  1  >  0 


d2  d  d 

(3.3)  mH -I(t)  +  cH(-^(t)-  ^.(t,())  +  kH(y(t)-u(t,  ())  =  0,  t>0 

dt 

respectively  where  y(t)  is  the  vertical  displacement  of  the  hose  mass  at  time  t  measured  from  the 
equilibrium  position,  f(t)  is  the  thruster  force  at  time  t  and  cH  and  kH  are  respectively  the  hose 
damping  and  stiffness  coefficients.  Assuming  that  the  rotational  inertia  due  to  the  hose  -  thruster 
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assembly  is  negligible,  we  also  have  the  zero  moment  condition 
^2  -\2 

(3.4)  CDI  —  ^(t,4>)  +  El  -4(t,  i)  =  0,  t  >  0 

9x2  3t  3x2 

at  the  free  end.  At  the  clamped  end,  we  have  the  usual  geometric  boundary  conditions  of  zero 
displacement, 

(3.5)  u(t,0)  =0,  t  >  0 

and  zero  slope, 

(3.6)  (t,0)  =  0,  t  >  0. 


Assuming  that  the  system  is  initially  at  rest,  we  have  the  temporal  boundary  conditions  or  initial 
conditions  given  by 


(3.7)  u(0,x)  =  0, 


du 

-^-(0,x)  =  0, 


0  <  x  <  l 


and 

(3.8)  y(0)  =  0,  $(0)  =  0. 

dt 


Our  primary  concern  here  is  the  inverse  problem  which  is  naturally  associated  with  the 
mathematical  model  given  by  equations  (3.1)  *  (3.8)  above.  The  physical  dimensions  and  mass 
properties  of  the  beam  and  thruster  assembly  and  the  elastic  properties  of  the  material  from  which  the 
beam  is  made  are  known.  Also,  the  thruster  output  can  be  experimentally  calibrated.  Consequently 
the  parameters  l ,  p,  E,  nvj.  and  b  and  h  (and  therefore  I)  are  known.  The  input  function  f  is  given  as 
well.  However,  the  coefficient  of  viscosity  cD  and  the  hose  parameters  mH,  cH  and  kH  must  be 
determined  via  an  identification  procedure.  Recalling  that  the  structure  is  instrumented  with  a  linear 
accelerometer  at  the  tip  of  the  beam,  we  formulate  the  following  inverse  or  parameter  identification 
problem. 

Given  a  known  input  f(t)  and  corresponding  measured  output  z(t)  for  t  e  [t^t,],  determine  the 
parameters  q  =  (cD,  mH,  cH,  kH)T  in  a  closed  and  bounded  subset  Q  of  R4,  which  minimize  the  least 
squares  performance  index 
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(3.9)  J(q)  =  f  I  ^-(U;q)-  z(t)  I2  dt 
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where  u(-,  •  ;  q)  denotes  the  solution  to  the  initial  boundary  value  problem  (3.1)  -  (3.8) 
corresponding  to  the  choice  of  parameters  q  =  (cD,mH,  cH,  kH)  e  Q. 

Implicit  in  the  statement  of  the  problem  above  is  the  well  posedness  of  the  initial  boundary  value 
problem,  i.e.  the  existence,  uniqueness  and  regularity  of  solutions  to  the  system  (3.1)  -  (3.8).  The 
infinite  dimensionality  of  the  distributed  state  constraints  necessitates  the  development  and  use  of 
some  form  of  finite  dimensional  approximation.  Both  of  these  issues  are  most  efficiently  addressed 
via  an  abstract,  functional  analytic  formulation  of  the  system  (3.1)  -  (3.8). 

Define  the  Hilbert  space  H  =  R2  x  L2(0,4)  endowed  with  the  usual  inner  product  and  let  V  = 
((£,Ti,cp)e  H:  <pe  H2(0,t ),  cp(0)  =  D(p(0)  =  0,  T|  =  <p(£)).  Define  the  coercive  bilinear  forms  c 
and  k  from  V  x  V  into  R  by  c(<p,y)  =  cH(£  -<p(£))(X.  -\j/(£))  +  cDI  <D2<p,  D2\p>and  k(tp,\p)  = 
kH(C  *  9(  ■£))$■  -  y(£))  +  El  <D2<p,  D2\(r>  for  <p  =  (£,  cp(£),  <p)  e  V  and  =  (X,  y( l ),  \p )  €  V, 
the  operator  M  e  £(H,H)  by  M(£,ri,(p)  =  (mH£,  myT),  p<p)  and  set  F(t)  =  (0,f(t),  0  )  e  H.  Then  the 
initial  boundary  value  problem  (3.1)  -  (3.8)  can  be  rewritten  in  weak  form  as 

(3. 10)  <  M  uu(t),  <p  >H  +  c(  u((t),  (p  )  +  k(  u(t),  $  )  =  <  F(t),  tp  >H,  t  >  0,  <p  e  V 

(3.11)  u(0)  =  0,  u,(0)  =  0, 


m 


for  u(t)  =  (y(t),  u(t.-fc),  u(t,  •))  e  V. 

Depending  upon  the  degree  of  smoothness  imposed  upon  the  input  f  as  a  function  of  t  (i.e.  Lj, 
Holder  continuity,  H1,  etc.),  standard  results  from  the  theory  of  abstract  parabolic  equations  (see  [F], 
[S])  can  be  used  to  demonstrate  the  existence  and  uniqueness  of  solutions  to  (3.10)  -  (3.11)  with 
varying  degrees  of  regularity. 

We  define  finite  dimensional  approximations  to  the  system  (3.10),  (3.1 1)  using  a  cubic  spline 

based  Galerkin  scheme.  ForeachN=l,2 .  let  denote  the  usual  cubic  polynomial  B- 

splines  defined  on  the  interval  (0, -6  ]  with  respect  to  the  uniform  mesh  {0,t/N,  2£/N,  ....  £  j  and 
which  have  been  modified  to  satisfy  pH(0)  =  DPj^O)  =  0,  j=l,2 . N+l.  Set  P|jf  =  (1,0,0),  PjN= 

(0,  Pj  (i ),  pj)  ,  j  =  1,2 . N  +  1  and  V  =  span  ( p^  .  The  Galerkin  equations  corresponding  to 

(3.10),  (3.1 1)  for  u^t)  g  VN  are  given  by 


(3.12) 


(3.13) 


<  Mu^(t),  pf>H  +  c  (^(t),  pf)  +  k(^(t),  P|*) 

=  <  F(t),  0H  >H,  t  >  0,  j  =  0,1,2 . N+l 


G"(0)  =  0 


u**(0)  =  0. 


WV ■  r 


For  q  =  (cD,  mH,  cH,  kH)  €  Q,  let  u^t;  q)  =  (yN(t;q),uN(t,£;q),  uN(t,-;  q))  denote  the  solution  to  the 
finite  dimensional  system  (3.12),  (3.13)  corresponding  to  the  choice  of  parameters  q  e  Q  and  let 

a2u 

u  (t;q)  be  defined  analogously  with  respect  to  the  system  (3.10),  (3.1 1).  If — -  (t,£;q)  in  the  least 

at2 

a2uN 

squares  performance  index  J  given  by  (3.9)  is  replaced  with  -  (t,£;  q)  and  the  resulting  finite 

at2 

dimensionally  constrained  estimation  problem  (henceforth  referred  to  as  the  N1*1  approximating 
estimation  problem)  is  solved,  an  approximation  to  the  optimal  parameters  q  is  obtained.  Indeed,  it 
can  be  argued  (using  estimates  in  the  spirit  of  those  given  in  the  previous  section;  see  [BGRW]  for 
the  details)  that  if  f  is  smooth  and  if  {qN}  ^=1  is  a  sequence  in  Q  with  qN  ->  q° , 

-.2  N  -.2 

then  u*(-  ;  qN)-»  uu(-  ;  q°)  in  L2(0,T;  H)  (or  equivalently,  that  — —  (•,£;qN)  -»  — -  (•,  l ;  q°) 

at2  at2 

in  L^(0,T))  as  N  -»  w.  This  in  turn  can  then  be  used  to  show  that  if  qN  is  a  solution  (which  can  be 
shown  to  exist)  to  the  Nth  approximating  estimation  problem,  the  sequence  {qN  }”=1  cQ  admits  a 

_N  <k>  ,  _Nj  _ 

convergent  subsequence  {q  }^=]  with  q  -4  q  as  j  — »  °°  and  q  is  a  solution  to  the  original  infinite 
dimensionally  constrained  identification  problem. 


Figure  3.4 
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We  tested  our  model  and  scheme  using  data  produced  by  the  following  experiment  carried  out 
for  us  on  the  RPL  structure  by  Dr.  Michel  A.  Floyd,  formerly  of  the  Draper  Laboratory.  With  the 
central  hub  held  stationary  and  the  structure  initially  at  rest,  a  thruster  line  was  pressurized  to  300  psi 
and  the  corresponding  jet  fired  for  50  milliseconds.  Linear  acceleration  at  the  tip  of  the  corresponding 
appendage  was  recorded  over  the  period  0  to  5  seconds  with  a  sampling  interval  of  5  milliseconds 
yielding  1000  measurements.  With  the  accelerometer  calibrated  to  5  volts/g  (g  =  32  ft/sec2),  the 
observations  produced  are  plotted  in  Figure  3.4  above. 

From  an  FFT  of  the  data  we  found  the  first  three  frequencies  to  be  approximately  .75  Hz,  7.5Hz 
and  MHz.  A  visual  inspection  of  the  data  immediately  reveals  the  modes  at  .75Hz  and  MHz.  The 
mode  at  MHz  is  a  torsional  mode,  or,  more  pecisely,  a  twisting  of  the  beam  about  its  longitudinal 
axis.  The  excitation  of  this  mode  results  from  a  combination  of  factors  including  the  presence  of  the 
hose  and  the  fact  that  it  is  connected  to  the  thruster  assembly  at  the  top  rather  than  the  center  of  the 
beam  and  the  opening  and  closing  of  the  thruster  valves.  Since  the  torsional  vibration  is  at  a  much 
higher  frequency  than  are  either  of  the  first  two  transverse  vibrational  modes,  and  since  it  is  relatively 
rapidly  damped,  in  our  study  here  we  simply  treated  it  as  noise  and  left  it  unmodeled.  A  more 
detailed  discussion  of  the  torsional  effects  and  its  coupling  into  the  accelerometer  measurements  of  the 
transverse  motion  of  the  beam  can  be  found  in  [BGRW]. 

The  beam  is  made  of  a  grade  of  aluminum  having  linear  mass  density  p  =  .027  slug/ft  and 
modulus  of  elasticity  E  =  15.84  x  108  lb/(ft)2.  We  have  l  =  4ft  and  I  =  bh3/12  =  4.71  x  10’8(ft)4. 
The  mass  of  the  thruster  assembly  nvj.  was  determined  to  be  .149  slug.  A  hose  pressure  of  300  psi 
was  determined  to  be  equivalent  to  a  force  of  .2971b.  We  have  therefore 

r 0.2971b  0  <  i  <  .05 

~  1 0.0  .05  <  t  S  5.0- 

The  index  of  approximation  N  was  taken  to  be  4  throughout.  A  detailed  study  of  the  convergence 
properties  of  our  scheme  using  test  examples  and  simulation  data  was  carried  out  with  the  results 
having  been  reported  in  [BGRW], 

We  neglected  the  hose  effects  and  internal  damping  (i.e.,  we  took  cD  =  mH  =  cH  =  kH  =  0)  and 
used  the  standard  Euler-Bernoulli  model  for  a  cantilevered  beam  with  tip  mass  to  generate  values  for 
the  linear  acceleration  at  the  tip.  The  resulting  acceleration  profile  along  with  the  experimentally 
observed  data  is  plotted  in  Figure  3.5.  If  one  is  willing  to  accept  the  Euler-Bernoulli  theory  as  a 
reasonable  description  for  the  transverse  vibration  of  a  long,  slender,  flexible  beam,  then  it  is  clear 
from  the  figure  that  the  hose  effects  are  indeed  significant  and  should  be  modeled.  The  residual  was 
found  to  be  3.03. 


Time  lSacar*3s) 


Figure  3.5 

In  applying  our  scheme  we  used  the  data  over  the  interval  3.0  to  5.0  seconds  (where  the 
contribution  from  the  torsional  mode  has  been  significantly  damped)  sampled  at  the  rate  of  10 
observations  per  second.  By  matching  the  first  two  modal  frequencies  of  the  data  with  the  first  two 
modal  frequencies  of  the  model  we  obtained  crude  estimates  for  mH  and  kH  that  could  be  used  as 
start-up  values  for  an  iterative  optimization  procedure.  Then,  using  our  spline-based  scheme  along 
with  an  iterative  Levenberg-Marquardt  nonlinear  least-squares  routine  to  solve  the  approximating 
estimation  problems  to  search  over  mH  and  kH,  optimal  values  for  the  hose  mass  and  hose  stiffness 
were  obtained.  Taking  these  values  along  with  cH  =  0  as  start-up  values,  we  used  our  scheme  again 
to  obtain  optimal  values  for  mH,  cH  and  kH.  Repeating  this  general  procedure,  we  obtained  the 
optimal  values  for  the  parameters  cD,  mH,  cH  and  kH  given  in  Table  3.1.  The  corresponding  tip 
acceleration  profile  is  plotted  in  Figure  3.6.  The  residual  was  computed  to  be  .70  -  a  clear 
improvement. 


cD  (lb  sec/ft2) 

mH  (slug) 

cH  fib.  sec/ft) 

kH(Wft) 

127.40 

.0801 

.0078 

.413 

Table  3.1 
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Figure  3.6 

In  summary,  our  results  demonstrate  the  feasibility  and  effectiveness  of  using  distributed 
parameter  systems  to  model  structural  dynamics.  Our  finite  dimensional  approximation  methods 
provide  a  viable  means  for  estimating  unknown  parameters  appearing  in  the  models  which  can  not  be 
explicitly  measured  experimentally.  Finally,  our  general  approach  appears  to  compare  favorably  (see 
[BGRW])  with  other  identification  procedures  for  the  estimation  of  parameters  in  models  for 
structural  dynamics  which  are  commonly  used  in  engineering  practice. 
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